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A mathematical analysis is presented for steady state heat conduction in cylinders, 
consisting of one or two isotropic materials disposed in concentric cylindrical volumes around 
the axis, heated by one or more continuous line heat sources parallel to the axis. The 
analysis considers both cylinders of infinite length, and cylinders of finite length subject 
to longitudinal heat flows arising from end conditions. Numerical solutions are given 
considering the line sources to be equal in magnitude and equally spaced angularly at a 
uniform radius from the axis. An outgrowth of this analysis has been the consideration 
of several new methods for the determination of the thermal properties of materials. 

1. Introduction 

[This paper presents a mathematical analysis of steady state heat conduction in a cylinder, 
consisting of one or two isotropic materials disposed in concentric cylindrical volumes around 
the axis, heated by one or more continuous line heat sources parallel to the axis. The analysis 
considers both cylinders of infinite length, and cylinders of finite length subject to longitudinal 
heat 1 flows arising from end conditions. The length of the line heat sources equals the 
length of the cylinder."] 

Problems concerning a single continuous line heat source are treated by Carslaw and 
Jaeger [l]. 1 The methods used in their book for single sources give an insight useful for the 
solution of problems involving multiple sources. In this paper, the multiple line sources are 
considered to be equal in magnitude, and to be equally spaced angularly at a specific radius 
from the axis. However, it is possible by the method used to obtain solutions for multiple 
sources having various and differing magnitudes, radii, and angular positions. 

This investigation was made in connection with an experimental steady state method [2] 
of determining the thermal conductivity of granular materials, using the arrangement shown 




Figure 1. Cross section of two concentric cylinders 
with line heat sources in the inner cylinder. 
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schematically in figure 1. The granular material comprises region 2. A long cylindrical 
ceramic core (region 1) is heated electrically by resistance wires in the equally spaced small 
longitudinal holes of the core. A thermocouple in the central hole of the core indicates its 
temperature at the axis, from which the average core surface temperature is inferred. If 
there is no heat source or sink inside the region enclosed by the circle of heater wires and 
there is negligible contact resistance between the two regions, it is intuitive that the tempera- 
ture in the hole of radius c at the center of the core is equal to the average temperature at the 
radius passing through the circle of heater wires. The holes are assumed to be virtual line 
sources. On this basis, a simple expression can be written for the average temperature at 
the radius of the hole at the center of the core : 
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where 

T c = temperature in hole of radius c at the center of the core, 
T 6 = uniform temperature at radius, r—b, 
9= heat input of line source per unit time per unit length, 
m= number of line sources, 
(r=hlk 2 , 

k x — thermal conductivity of inner cylinder material, 
k 2 = thermal conductivity of outer cylinder material, 
a = outside radius of inner cylinder, 
b= outside radius of outer cylinder, 
r'=radius passing through heat sources (r'<^a). 

If the other fixed quantities in it are known, eq (1) serves for the experimental determina- 
tion of a, or, if an adequate value <j can be assumed, of k 2 , the thermal conductivity of the 
material of the outer cylinder. 

In this paper, eq (1) is derived analytically for cylinders of infinite length, and an exam- 
ination is made of departures if the core temperature were measured not at the center of the 
core, but at a small radius and variable angular position. For cylinders of finite length, with 
temperatures measured at midlength, departures from eq (1) due to longitudinal heat flow 
caused by end effects are given for variations in a and cylinder length. 

2. Mathematical Development 

2.1. Infinitely Long Cylinder 

a. Single Line Heat Source 

Consider a single continuous line heat source of magnitude Q heat units per unit time per 
unit length and infinite in length, located in an isotropic medium at radius r' , polar angle cp' . 
At a point in the medium, (r, <p), the distance measured from the source is R=[r 2j rr' 2 — 2rr' cos 
(<p— <p')]K The heat flux from the source at a line of unit length parallel to the source must be 
satisfied by the limit 



lim 



[-^asj-fc (2) 



where k is the thermal conductivity of the medium containing the source and 6 is the tempera- 
ture above an arbitrary datum temperature. The above limit is satisfied if 

8=-£- k lnB. (3) 
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The following identities are necessary for further development of the problem 

i /> i / ^ /r\ n cos n(<p—<p') ,, 

=lnr-sCT c08 ^-^ r>r'. (4) 

71 = 1 \F/ ^ 

By use of the above identities, eq (3) can be readily shown to be a particular solution to a 
partial differential equation for heat conduction: 
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Here, the source is assumed to be contained in region s=l, where c<Cr<Cji; for region 6* = 2, 
a<jr<Cb. Solutions to (5) are: 

For region 1, c<V<V 

»■= A [-- ™+= '-^F^ {(?)"+-•- (£)+*■ (£)"}] 

For region 1, r'<Cr<^a 
For region 2 

*-"sk [" ln ' +i> '+S 25i»fa^ {C^+ZV-}} 

The constants Pi, P 2 , A, -4», #n, C n , and Z> n are determined from the boundary conditions: 

2 =O at r=b 

^=0 atr=c 

or 

#i=0 2 at r=a 

k 1 ^=k 2 v^ at r=a. (7) 

The boundary condition, ^r = at r==c assumes that the region 0<r<c contains a perfect 

insulator, and does not assume the temperature at r=c to be constant. This boundary condi- 
tion is a mathematical expediency for the development to the problem of more than one line 
heat source, where the heat flows directed toward the center from heat sources positioned on r' 
tend to nullify each other at r=c. The degree to which the gradient becomes negligible at 
r = c depends upon the uniformity of the magnitudes and positions of the heat sources. From 
a practical viewpoint, this boundary condition approximates a system where region 1 is com- 
posed of a material of high thermal conductivity and the region 0<r<c is a thermal insulator. 
Otherwise, for the single heat source, the region 0<r<c, and its thermal conductivity must be 
introduced into the analysis. 

Substitution of the boundary and continuity conditions in eq (6) gives for region 1 

*-£ !.»«//+,>„ wH^gP-SfaQ^) •+(£)"+«-•)] <s) 
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for c</<V, where 

*.(r) = 



(*xz+i)[i+$ n f + W + W ^ J 



l-(a/6) 2 * 

l + (a/6) 
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For r'<Cr < Ca, interchange r and r'. 

b. Multiple Line Heat Sources 

Consider m continuous line sources of magnitude Q p at (r' } <p' p )j p=l, 2, . . . , m. The 
separate solutions involving <p' p are additive. Thus, 

! s ! - (p)"+G£) S +*»W - 
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If Q P = Q, and <p p = 2Tp/m, i.e., equally spaced equal line sources, then 
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=0 for q^mn 

=m cos m^ for q=mn, 
and eq (10) becomes 2 

where q=mn. 

Using (4), and rearranging, eq (11) becomes valid for the whole region 1, c<V<a, 

-| In |l+(^) 2W -2 (£)" cos mv |+ mS ^|*?^(r)]. (12) 

Inspection of (9) shows that the ratios X, Z, and the dimension ratios are less than unity, 
so that the expression \p g (r), in (12), for m>2, becomes quite small with increase in n, and 
only a few terms of the series need be evaluated. 

The average temperature, 6 T , at any radius r, can be found by integrating eq (11) with 
respect to <p, or 

e r =^ jr J e x d<p=^ [In a/r+a In b/a], r'<r<a 
= £f [In a/r' + v In 6/a], c<r<r' . (13) 



2 Equation (11) may also be derived by considering a wedge, t=re i ?, 0O<27r/7ra, c<r<b, with a line source at r=r' lying on its bisector. 
Due to the symmetry of the line sources assume no heat flow across the boundaries of the wedge. Making the transformation w=t m gives a single 
line source at r = r' m and <p'=ir. Using the single line source solution (8), a transformation from the w to t plane gives (11). 
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Equation (13) substantiates eq (1), which was arrived at intuitively. Similar integration of 
(10) with respect to cp leaves only the first term of the right side of (10). If it is assumed 

m 

thai mQ a v=S Qpj then eq (1) is valid for reasonable variations of the distribution and mag- 
p-i 

nitude of the heat sources on the radius r=r r (see discussion following (7)). 

For a finite number of sources, it is of interest to investigate the ratio 

fiJf (14) 

at r=c, as a function of the angle <p for variations in the number of line sources and in a, the 
ratio of thermal conductivity in the first region to that in the second. Numerical solutions are 
shown in figures 2 and 3, where the radial dimensions were fixed by means of the following 
ratios: 6/o = 5.0, r'/a = 0.7, and c/a=0.2. 
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Figure 2. Infinite composite cylinder — Ratio eq 
(14), versus angle for increasing number, m, of line 
sources, with a= 10, b/a = 5, r'/a = 0.7, and c/& = 0.2. 
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Figure 3. Infinite composite cylinder — Ratio (3 eq 
(14), versus angle for various values of <r( = ki/k2), 
with m=3, b/a = 5, r'/a=0.7, and c/a = 0.2. 



2.2. Cylinder of Finite Length 

a. Single Line Heat Source 

The heat flux from a continuous line heat source in the finite cylinder (0<2<7)7must 
satisfy the limit (2). A solution that satisfies (2) and the partial differential equation 



^LlI ^_j_I ^4_^ = n 



(15) 



is 



20 - sin T~ 

TT 2 fC w = i J 



<¥) 



j=2n-1. 



(16) 
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The following identities and definitions are necessary for further development of the problem 
Ko(«#)= £ cos ifa-rfKiiar'Vtar) r<r' 
K (aR)= £ cos Hv-rtKAcrMar') r>r' 

i= — co 

F„(x, y)=I n (x)K n (y)-I n (y)K n (x) 
F' n (x,y)=I n (x)K' n (y)-I' n (y)K n (x) 

F';{x, y)=r n {x)K' n (y)-r n {y)K' n {x) 
F n {x,x)=F';(x,x)=0 
F' n (x, x)=-l/x. (17) 

Where I„(x) and K„(x) are modified Bessel functions of order n, and of the first and second 
kind, respectively. 

Proceeding as in section 2.1. a, and using (17), the solutions for the two finite regions of 
the composite are: 



for region 2, a<Cr<^b 

2Q ^° 
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for region 1, r'<V<^ 
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cos «(^V) {^ (f ) /, (^)+^ (*? )} 
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for region 1, c<><V 



^=H fj ^ .£ cos i (,-,') { Ki (*f)l t (if) +Ui (ff)} 



where 



*(^a (*?>■«*(*?} 



M z -, Pj, and Q* must satisfy the boundary and continuity conditions (7). The exterior bound- 
ary condition, that 0=0 at the surfaces 2=0, z=l, and r=b, is satisfied by the above equations. 
For region 1, c<V<r' 

20 Sm V 

^=zir S — t— S y«( r 'i r ) cos ^-V) (19) 



where 



i- ; (if, g) [> ; (M, g) ,, (y , a)_^ ; (y, g) , ; (y, gy| 

7,(r , r) = - 

w* (M, if) F , (t» i-yr, (if, if) F[ («, fe) 
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b. Multiple Line Heat Sources 



Consider, as in section 2.1.b, m equally spaced equal sources at r=r'. Equation (19) 
becomes 



Or 



. mz 
V S — r- < 7o(r',r)+2 ^ T.fr',r) cos qy )- 



% 



U(r',r)+2g 7 ,(r',') cos 2v J- 



(20) 



where q=mi andi=2n-l. 

It is of particular interest to evaluate the mean temperature at r=c and 2=1/2, 

- 2mQ - (-1)» +1 , , , 
(0c)z=i/2=—2T S 7 7o(r ,c) 

to determine the effects of longitudinal heat flow by the ratio, r, of this temperature to the 
average temperature for the infinitely long cylinder at r=c as given by (13): 



(o t ) 



z=U2 



".■ 



(21) 



Numerical solutions for r versus the dimension ratio l/b for various values of a, are shown in 
figure 4 for fixed radius dimension ratios b/a=5.0, r f /a=0.7, c/a=0.2. 

I.Or 



Figure 4. Fin /'(<■ composite cylinder — Ratio t eq (21) 
versus l/b (cylinder length/outer radius) for various 
values of a, with b/a = 5, r7a = 0.7, and c/a=0.2. 




2.3. Special Cases With cr=l 

An interesting special case of multiple line sources in an infinitely long cylinder is the 
case where k\=k 2 . Substitution of a=l in (11) gives 






^^m<^'-^>] 



(22) 



where 



#M= 



'+C);(^);(£) 



For ?'=c and for c<<r'<6, only the first term in (22) need be considered. Physically, this 
special case represents a single cylinder of homogeneous material, in which m heater wires 
parallel to the axis are embedded at a radius r' with equal angular spacing. The temperature 
is measured at the axis of the cylinder. 

For the special case in a finite cylinder where c = and k\ = k 2 or <r=l, the y t of eq (19) 
becomes 



y t (r', r)= 



T (kirr\ T1 /kwb kirr'\ 



'•(") 
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3. Discussion 

The variations in temperature at a radius c (c<V) with respect to angle are illustrated 
in figures 2 and 3, using fixed values of the radius dimension ratios b/a, r'\a and c/a. Figures 2 
and 3, for the infinitely long cylinder, show that the angular variation decreases with increase 
in the number of line sources (?n), and also with increase in the value of a= (k\lk 2 )- The latter 
is intuitively evident since the magnitudes of all temperature differences in the inner cylinder 
obviously are reduced as its conductivity increases relative to that of the outer cylinder. 

Figure 4 shows the above-datum temperature at the center and midpoint of a cylinder 
of finite length as a fraction of the corresponding temperature attained in an infinitely long 
cylinder (eq (1)), for variations of a and of l/b, the ratio of the cylinder length to the outside 
radius of the outer cylinder. Since this analysis assumes that the ends of the two concentric 
cylinders are maintained at the same temperature (6=0) as the outer surface at r=b, the l/b 
ratio necessary to maintain a specified r-value from figure 4 may be somewhat different from 
that encountered in most applications. With an increasing value of a, longitudinal heat flow 
in the inner cylinder increases relative to the radial heat flow, hence it is necessary to increase 
the cylinder length to avoid an increase in the error. For values of o much less than unity, 
the length of cylinder necessary to restrict the error is controlled chiefly by the effect of lon- 
gitudinal heat flow to the ends of the outer cylinder. 

This analysis confirms that for steady radial heat flow in angularly isotropic infinitely 
long circular cylinders, the temperature at the axis is equal to the average temperature over 
any coaxial cylindrical surface which does not envelop a heat source. Thus, in such a cylin- 
drical system having steady line heat sources parallel to the axis and all at a radius r', the 
temperature at the axis is equal to the average temperature at the cylindrical surface of radius r' , 
If the line sources were uniformly spaced angularly, there would be no difference in average 
temperature or heat flux at a radius greater than r f , whether the heat was supplied by the 
distributed line sources and the average temperature at r' was determined by a measurement 
at the axis, or whether the same aggregate heat was supplied at the axis, and the temperature 
at r' was determined by averaging the temperatures measured at the several angular positions 
at r'. In short, the axial position and the equiangular positions at r' ', respectively, can be 
used interchangeably as sites for heat supply or temperature measurement, as convenient, 
without affecting the average heat flow pattern in the cylinder beyond r' . 

The use of equiangularJy distributed line sources for heat supply has advantages in some 
applications. For example, the apparatus described by Flynn [2] used a cylindrical ceramic 
core as a central heat source in testing a hollow cylinder of fine granular material. If the 
granular material were very coarse, the thermal contact between the granular particles |and 
the smooth convex ceramic surface might be enough different from that between particles to 
cause an unwanted and not easily ascertainable temperature drop in the critical region near 
the ceramic core. If, on the other hand, such a coarse granular material were tested as a 
complete cylinder, with temperature measured at the axis and line heat sources at radius r' in 
the material (i.e., the case of <j-=l, eq (22)), the oddities of thermal contact of particles with 
the line heat source would be inappreciable at short distances from the line sources, since the 
heat flow there would be that occurring in the particulate region. 

Another example, of possible experimental interest, is that of conducting a measurement 
of the thermal conductivity of a solid by the Powell "stacked-disk" method [3], but using 
the axial hole for the temperature measurement and the inner circle of holes (preferably three 
or more) for heat supply, in reverse of the arrangement used by Powell and by others. In 
this way, uncertainties arising from measurements of temperature in regions of large tempera- 
ture gradient, or due to disturbances of the heat flow pattern by the thermocouples, would 
be largely avoided. 



126 



The author is indebted to J. C. Jaeger, Professor of Geophysics in the Australian National 
University, for his comments on an earlier draft of this paper. 
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